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Abstract 

The occurrence of a big smash singularity which ends the universe in a finite 
time in the future is investigated in the context of superquintessence, i.e. dark 
energy with effective equation of state parameter w < — 1 and H > 0. The 
simplest toy model of superquintessence based on a single nonminimally coupled 
scalar field exhibits big smash solutions which are attractors in phase space. 
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A new picture of the universe has emerged in recent years: the data from the Boome- 
rang [1] and MAXIMA [2] experiments confirm that we live in a spatially flat (Q = 1) 
universe, where f2 = is the total energy density expressed in units of the criti- 

cal density [3]. Baryonic and dark matter only account for fi( m ) ~ 0.3 of the total energy 
density Q, while the rest is due to a yet unknown form of dark energy. Studies of Type la 
supernovae [4] and of radio-galaxies [5] show that the present expansion of the universe is 
accelerated, i.e. a > 0, where a(t) is the scale factor of the Friedmann-Robertson- Walker 
line element 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ^) (1) 
describing our universe in comoving coordinates (t,x,y,z). The Friedmann equation 

- = -- p + 3P , (2) 
a o 

where p and P are, respectively, the energy density and pressure of the source of gravity, 
shows that in order to have accelerated expansion the pressure of the dark energy do- 
minating the dynamics must be negative, P < — p/3. A cosmological constant A as the 
explanation of dark energy is rejected by most cosmologists because of the cosmological 
constant problem [6] and of the cosmic coincidence problem. The vacuum energy density 
predicted by high energy physics with a Planck scale cutoff is wrong by 120 orders of 
magnitude (or 40 orders of magnitude if the cutoff is at the QCD scale). Solving the 
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cosmic coincidence problem of why the dark energy came to dominate the dynamics only 
recently (at redshift z ~ 1) requires extreme fine tuning of A. One would rather have A 
exactly equal to zero due to some yet unknown mechanism than this extreme fine-tuning, 
and it is preferable to explain the present cosmic acceleration with a dynamical vacuum 
energy called quintessence. Many models of quintessence, most of which based on scalar 
fields, have been proposed. 

Observational efforts aim at determining the effective equation of state parameter of 
the universe w = P/p. It was pointed out that the current data allow for [7, 5, 8, 9], 
or even favour [10]- [15], values of this parameter in the w < — 1 range. Were a value 
w < — 1 to be confirmed by observations, it would be a most interesting finding, because 
such values cannot be explained by Einstein gravity with a minimally coupled scalar, 
if one assumes positivity of the energy density. In fact in such "canonical" models the 
scalar field has energy density and pressure 

P=^- + V(4>), P= L f-V(t), (3) 
and effective equation of state parameter 

w = —i- (4) 

where x = (J>) / (2V) > if V > 0, giving — 1 < w < 1 (the minimum of w being 
attained by de Sitter solutions). In these models one only obtains w < — 1 by assuming 



V < — (0J /2 < 0, which in turn implies a negative energy density (3), and p < is 
hardly an acceptable proposition. Correspondingly, the Friedmann equation 

H = -^(p + P) (5) 

tells us that, for a minimally coupled scalar in Einstein's gravity it is H = —k (J^j /2 < 
(the extreme case H = again corresponding to de Sitter solutions). A regime with 
w < — 1 is associated to H > and is called superacceleration; a form of dark energy 
capable of sustaining superacceleration was dubbed superquintessence [16] or phantom 
energy [10]. 

Models have been proposed to explain superacceleration regimes, including scalar 
fields nonminimally coupled to the Ricci curvature, actions with the "wrong" sign of 
the kinetic energy of the scalar, supergravity-inspired models with non-canonical kinetic 
energy terms and zero potential (k-essence), Brans-Dicke-like fields in scalar-tensor gra- 
vity, or stringy matter [17]. If the universe superaccelerates, its expansion becomes so 
fast that it risks ending its existence in a finite time, a(t) — > oo as t — > to with to finite. 
For solutions with this property (big smash or big rip solutions) [10, 18, 19, 12, 14] the 
energy density of superquintessence increases with time instead of redshifting away as 
the matter or radiation energy densities oc a -3 , p^ oc a~ 4 , or as the energy density 
of "ordinary" quintessence. Big smash solutions occur in the theory of ordinary diffe- 
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rential equations (ODEs) when the solutions of an ODE cannot be maximally extended 
to an infinite interval. An example of ODE which does not admit maximal extension of 
its solutions is 

— — = Ay (x) ; (6) 

if A > the rate of change of the solution is fed by the increasing value of y itself in 
a positive feedback mechanism that makes y(x) grow so fast that it explodes in a finite 
time, while this behaviour is absent if A < 0. As we shall see, eq. (6) is similar to the 
equation satisfied by the Hubble parameter of a superaccelerating universe. 

A big smash can be avoided in certain models of superquintessence (e.g. [11]) but 
is a generic feature in other models [10]. Whether the big smash is unavoidable or not 
depends, of course, on the model adopted [20]: here we consider what is perhaps the 
simplest model of superquintessence, namely a single scalar field coupled nonminimally 
to the Ricci curvature, described by the action 

S = J d'x^g [| (i - £0 2 ) - \<f* V c <f> V d <j> - V{4>)] + , (7) 

where £ is a dimensionless coupling constant and is the action for ordinary matter. 
Nonminimal coupling is introduced by renormalization even if it is absent at the classical 
level and is required in classical general relativity by the Einstein equivalence principle 
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[21]. The field equations are 



G ab = k T ah [(f)] , 



(8) 



T ab [0] = V a V 6 - - <uV c V c - V g ab + £ - V a V 6 ) (0 2 ) + £G afe 2 . (9) 

The gravitational coupling in eq. (8) is the usual constant k = 8nG and not the ef- 
fective time-dependent coupling K e ff — k (1 — k£0 2 ) 1 recurring in the literature and 
corresponding to a different way of writing the field equations [22]. Moreover, the scalar 
field stress-energy tensor (9) ("improved energy- momentum tensor" [23]) is covariantly 
conserved, V b T ab [4>\ = 0. In the metric (1), the field equations become 
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dV 



1 - f (1 - 60 /t0 2 l (H + 2H 2 ) - k (6f - 1) 2 - AkV + 6<0 — = , (10) 

d(p 



2 + 6£ K,H(f)(j) - 3H 2 (l - <0 2 ) + kV = , 

+ 3if + £i?0 H — 7— = . 

a0 



(11) 
(12) 



Only two of these equations are independent, as the Klein-Gordon equation (12) can be 
derived from the conservation equation + 3H ^p(^) + = when ^ 0, where 



P W = y + ^(0) + 3^0 (H<P + 20) 



(13) 



and 



2 



P w = y - ^(0) - f [4/700 + 20 2 + 200 + (2if + 3H 2 ) s 



(14) 



are, respectively, the effective energy density and pressure of the fluid equivalent to the 
nonminimally coupled scalar. Note that the Hamiltonian constraint (f f) can be written 



as 

H 2 = ^P^, (15) 

which is consistent with p^ > 0. This is not the case with other definitions of effective 
energy-momentum tensors T afe [0] used in the literature (see the discussion in Ref. [22]). 
In order to investigate the fate of the universe with respect to big smash solutions, 
we neglect the matter part of the action S^; this assumption is justified by the fact 
that, when a superacceleration regime sets in, the superquintessence energy density 
quickly grows to dominate the matter energy density, which instead fades away. 

Only the two variables H and are needed to describe the dynamics of the system 
(10)-(12), and the phase space is a two-dimensional manifold with a rather complex 
structure [24]. This is best seen by rewriting the field equations as [24] 



W±-^(F^), (16) 



H 



3 (2f - 1) H 2 + 3£ (6f - 1) (4£ - 1) kH 2 4> 2 + £ (6£ - 1) v^#0 

1 



dV 

+ (1-20 kV-kS<I> — 



l + <(6£-l)0 2 ' 



(17) 
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where 



T (H, 0) = 8k 



2 \3H 2 



V(</>) + 3£ (6£ - 1) H 2 4> 2 



(18) 



K 



The appearance of the ± signs in eqs. (16) and (17) requires a clarification: the phase 
space curved manifold is composed of two sheets, corresponding to upper and lower sign, 
and there typically is a region forbidden to the dynamics corresponding to T (H, <p) < 0; 
the two sheets join each other at the boundary of this forbidden region, which needs 
not be simply connected. Orbits of solutions lying in the "upper" sheet can switch to 
the "lower" sheet at these points, and vice- versa [24]. Far from the forbidden region's 
boundary, the orbit of a solution is forced to stay in one sheet and cannot cross over 
to the other sheet without approaching the forbidden region and touching its boundary. 
Hence, for large values of H and only one sign in eqs. (16) and (17) applies to each 
solution. Since H(t) and (f>(t) grow so quickly when the superacceleration regime sets 
in, it is meaningful to perform an asymptotic analysis for large values of these variables 
when searching for big smash solutions. We consider conformal coupling £ = 1/6, which 
is a stable infrared fixed point of the renormalization group [25], and consider as a toy 
model the potential 



where A will be required to be negative. This does not harm the positivity of the energy 



v{4>) 




+ A0 4 , 



(19) 
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density when £ 7^ 0. Eqs. (16) and (17) then reduce to 



0~ -#0± V^2\<i) 2 , (20) 

if ~ -2H 2 + fMj) 2 , (21) 
where \i = Km 2 /6. We look for big smash solutions of the form 

a(t) = j^h^ > ( 22 ) 



|t-to| c 



and 



= u Vife > (23) 
I ^ — *o I 

with a±, f3± > 0, consistently with the approximation of large H and (f> employed, and 
where to, and 0* are constants; to will be approached from below. Substitution of 
eqs. (22) and (23) into eqs. (20) and (21) yields 



±y/-X (2/x + A) —(fX + X) 

a ± = VTlx ■ (24) 

f3 ± = 1 , (25) 

By taking the positive sign in eq. (24), one immediately sees that there are big smash 
solutions (a > 0) in the range of parameters 

K ttt < 4 • (27) 



Next, one would like to know whether these big smash solutions are stable or if they 
disappear when perturbed. The equations for the perturbations 5H and Sep are suffi- 
ciently involved to defy a direct analytical investigation. It is convenient, instead, to 
focus on the projection onto the (if, 0) plane of the two sheets composing the phase 
space. Assuming H, <fi ^ 0, as is legitimate during the late stages of a superacceleration 
regime, one obtains the vector field 

dH H 2u 2 - /I 



d(f) u =F V-2A 

where u = Hj 0. The identity 



(28) 



dH , du ,_ x 

#=" + *# (29) 



then yields 



du = u 2 ± V-2\u- fi ( l\ 
d<P ~~ u Tv QI '{(/>)' 

It is straightforward to see that exact solutions of eq. (30) with u =const.^ ±V— 2A 
exist and are given by 

Eq. (31) includes the big smash solutions (22)-(26). In fact, simple algebra shows that 
for the latter 

a ± y/4/i - 2A + y^2A 

u= vr 2 ■ (32) 
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which reproduces a special case of eq. (31). 

Let us proceed to study the stability with respect to linear perturbations c 
solutions of eq. (30), which can be rewritten as 

H = 7i , 

where the constant 7$ can assume the values 



-x/^2A + V4/i-2A ^ 

71 = 2 >0 ' 

-v 73 ^ - V4/i - 2A n 

72 = g <0, 

73 = -71 < , 

74 = -72 > 

(the big smash solutions (22)-(26) corresponding to 72). The perturbations 8u in 

= m + Su = 7j + 5-u (0) , 

satisfy the equation 



d (<5w) 2-u 5m ± \J — 2A 5m 
-u =p V-2A 



which yields 



where e is a constant and 



Su = e (f) 1 , 



. 2 7i ± v^2A 
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There are the following possibilities for 5i, corresponding to the four values of 7$ and to 
the upper and lower sign in eq. (41): 

which is positive if fi/ |A| > 17/12 ~ 1.4167; this range of values of the parameters 
corresponds to growing perturbations and instability, while the finite range < /// |A| < 
17/12 corresponds to stability. 



f u = 2 ^-!±S (43) 
is positive if ji/ |A| > 7/24 ~ 0.2917; this range corresponds to instability. 



_ 2^4^ - A 

^2a = ; — ; < (44) 

corresponds to the big smash solutions and to stability for any value of the parameters 
A < and \x\ on the other hand, the possibility 



V /4^A + 2 V /^2A 

52b = 2 7^W^T (45) 

corresponds to instability. The remaining cases give the values of 5 



^3a — 8 lb , 5 3b — 8i a , 

ha = 5 2 b , S 4 b = 5 2a (46) 
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already considered. 

The big smash solutions are stable against linear perturbations and behave as at- 
tractors in phase space. Thus, there is a finite chance that a universe described by the 
model considered here end its existence in a finite time due to a big smash with infinite 
expansion, in which the energy density diverges instead of being diluted away and bound 
systems are gradually ripped apart [14]. 

A few considerations are in order: first, the value of the effective equation of state 
parameter w is still subject to uncertainty and a value w < — 1 associated to superac- 
celeration is not yet confirmed; second, even if such a value were supported by future 
experiments, it does not automatically imply that the universe will end in a big smash. 
Third, the toy model of superquintessence employed here should be generalized to more 
realistic models, which however are not constrained sufficiently well by the presently 
available data. Potentials used to model quintessence are usually different from (19), 
but the latter is a very common potential in scalar field cosmology and is easier to study 
as a toy model (it is difficult to reach conclusions about the existence of big smash so- 
lutions with other potentials, and even more difficult to perform a stability analysis). 
The presence and stability of big smash solutions in a finite future for a wide range of 
parameters leads one to regard a big smash as a generic feature of scalar field models 
of superquintessence that include a nonminimal coupling to gravity. This can be of the 
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simple form described by the action (7), or of a more general form as in scalar-tensor the- 
ories, which have been known to contain smash solutions for a long time [26]. The action 
(7) can be explicitly reformulated as a scalar-tensor theory with variable Brans-Dicke 
parameter 

^ = mr^i • (47 > 

and 

V = q , (48) 

but more general forms of the coupling function u(4>) are possible. Since the field 
equations for the coupled variables H and (f> in scalar-tensor gravity exhibit terms similar 
to the right hand side of eq. (6), solutions with explosive growth are possible in these 
theories. 

The most stringent constraint on the theory of the nonminimally coupled scalar comes 
from Solar System experiments. The Brans-Dicke-like field ip mediates a long range force 
that is constrained by tests of general relativity. Since <p varies on a cosmological time 
scale, it is appropriate to approximate (p with its present value <p and uj (ip) ~ uj (<p ) = 
lu . The lower bound on uj is 1 \uj \ > 500 [29, 30], which yields 



0J 



2 



> 500 . (49) 



1 The more stringent bound \uo\ > 3300 [31] yields a constraint on |^| of the same order of magnitude. 
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Although the present day value of is unknown, a weak coupling regime in which 
k |£| 0q << 1 is plausible [27], given that typical values of £ predicted by renormalization 
are of the order of 10 _1 -10 -2 , and that O cannot exceed the Planck mass m p i = G -1 / 2 by 
too much without causing fine-tuning problems in the parameters m and A (the energy 
scale V must be below the Planck energy scale m pl ). In the weak coupling regime 
k|£|0q << 1 and assuming that y/n |0 O | — 1 (corresponding to O — 0.2m pi ), one 
obtains |£| < 5 ■ 1CT 4 (this limit is weakened if \<f>\ « m p i). This constraint limits the 
amount of superacceleration that is present today (see e.g. Ref. [28]), but it should be 
kept in mind that even a small amount of superacceleration can eventually lead to a big 
smash. Stringent limits on |£| can be satisfied and the big smash will occur later: whether 
this amounts to fine-tuning the value of £ is determined by the still unknown value of 
the parameter w which quantifies the present amount of superacceleration (assuming 
that the universe really does superaccelerate, i.e. that w < — 1). It is hoped that the 
observational determination of the value of w will soon clarify this issue. 

To conclude, even if the departures of gravity from general relativity are small to- 
day in the Solar System, on a large scale they may have a catastrophic effect on the 
future of the universe. Usually, research on scalar-tensor cosmology has focussed on 
how scalar-tensor theories can depart from Einstein's gravity in the early universe and 
converge to it at later epochs. Here, the issue is rather the one of a cosmological solu- 
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tion of scalar-tensor gravity that is close to a general relativistic solution today (when 
superacceleration is still moderate and w close to —1), but will dramatically depart from 
it in the future. If the universe really superaccelerates, the concern about a big smash 
is legitimate. It is intriguing that observational data place our present universe so close 
to the boundary w — — 1 between the possibility of a big smash and certain evolution 
into infinite dilution in an infinite time. 

The author acknowledges Leon Brenig for a stimulating discussion at a meeting in Peyresq. 
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